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Recent work has shown that boundary correlators in AdS-Schwarzschild 
can probe the geometry near the singularity. In this paper we aim to analyze 
the specific signatures of the singularity, show how significant they can be, 
and uncover the origins of these large effects in explicitly outside the horizon 
descriptions. 

We add perturbations to the metric localized near the singularity and 
explore their effects on the boundary correlators. Then we use analyticity 
arguments to show how this information arises from the Euclidean path in- 
tegral of a free scalar field in the bulk. 
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Introduction 

In physics we are often able to understand a new object by probing it with 
something we are already familiar with. This standard approach fails when 
applied to spacelike singularities because our probes cannot get close to 
spacelike singularities and live to tell about them - they become trapped 
behind horizons longs before they approach the singularity itself. However, 
the structure of a singularity can be probed by spacelike geodesies, and in the 
AdS/CFT correspondence such geodesies approximate opposite sided bound- 
ary correlators. These ideas helped motivate previous work pQ, [3], [3], and 
[5] where it was shown that information about the singularity could be ex- 
tracted from these correlation functions. Our purpose here is to show that 
these correlators really do contain information about the geometry near the 
singularity, and to give explicit examples demonstrating that it is 'easy' to 
access this data. Most of these examples utilize the geodesic approximation. 

In jl] it was shown that there exists a dual description of physics in 
D = 3 AdS-Schwarzschild, depending on whether one analytically continues 
the Euclidean metric with respect to Schwarzschild or Kruskal time. In the 
Schwarzschild description, all quantities can be determined with outside the 
horizon calculations. Since we are extracting information about the neigh- 
borhood of the singularity, it is surprising that we can see large effects from 
outside the horizon computations. Thus our final calculation will be an at- 
tempt to show how the large effects arise from the full quantum path-integral. 
The key will be analyticity. 

The layout of the paper is as follows. After a short review of the AdS- 
Schwarzschild spacetimes, we use the WKB approximation to explore the 
effects of a simple perturbation localized near the black hole singularity in 
the D = 3 case. This is simply a small change in the geometry near the 
singularity; it is not required to satisfy Einstein's equation and we do not 
expect it to be a model for anything physical. Its purpose is to illustrate 
that we can explore the geometry with boundary correlators by studying 
them as a function of the boundary time. We find that the perturbation 
produces a significant change in the correlators for relatively small values of 
the boundary time, indicating that information from behind the horizon is 
not difficult to extract, at least for a meta-observer who can compute any 
gauge theory amplitude. Moving on to a D = 5 example, we make a very 
naive 'stringy perturbation' to the metric that resolves the singularity, and 
investigate the consequences (this is not intended as a guess at what string 
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theory does, it is just a possibility). We find that a pole in the opposite sided 
correlators resulting from geodesies that 'bounce' off the singularity in the 
unperturbed case splits into two poles above and below the real axis. This 
is a distinct, unambiguous signature that we have removed the singularity. 

Previous work has focused on the idea of behind the horizon data made 
available from amplitudes computed outside the horizon, far from the black 
hole. Thus as a concrete example, we calculate approximately how many 
terms in a power series of the opposite sided correlator as a function of 
boundary time one must compute in order to see a large change resulting 
from a small perturbation at the singularity. We find that for a perturbation 
decaying as e~ ar it is only necessary to calculate 0(a) terms. Finally, we set 
up the full quantum calculation of correlation functions for a free scalar field 
in a perturbed Euclidean D = 3 background. We expand to first order in 
the perturbation and derive a geodesic approximation, showing that in the 
Euclidean case the quantum computation is consistent with the WKB ap- 
proximation. This shows by virtue of the analyticity of multi-point functions 
that the large effects of the perturbation must be reproduced in path integral 
computations. 



Spacetime Geometry 

Here we will be considering the effects of perturbations to the D = 3 and 
D = 5 AdS-Schwarzschild geometry localized near the black hole singularity. 
In both cases the metric takes the form 

ds * = _ f{r)dt 2 + dr^_ + rW (i) 

Henceforth, we will ignore spherical part. The function 

f(r) =r 2 -l- g{r) (2) 

for the D = 3 case and 

/(r) = r 2 - 1 - £(r) (3) 

for the D = 5 case, with g(r) a 'small' perturbation function for both cases. 
It is essential that we always choose an analytic g(r) function to maintain 
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Figure 1 : This figure shows the Penrose diagram for the D = 3 AdS-Schwarzschild 
spacetime with a symmetric spacelike geodesic passing near the singularity. The 
D = 5 case has a similar diagram, with one important difference - the singularities 
must 'bend outward', as shown in jSJ. 



equivalence with the boundary gauge theory, whose correlation functions are 
analytic. In these backgrounds the geodesic equations take the form 

E 2 + f(r) (4) 
E (5) 

The WKB approximation to the opposite-sided boundary two-point func- 
tion is given by e~ mL , where L is the proper length of the geodesic connecting 
the two points. Thus these geodesic equations will be used extensively, be- 
ginning in the next section. 



dt 
as 
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Some D = 3 Examples 



Geodesies Far From the Singularity 

First we will give a very simple argument showing in general that the effect 
of any g(r) perturbation is appropriately small for geodesies far from the 
singularity. Integrating equation (4) we find 



lim 

R c — >oc 



Rc dr 
2 / — = . - 2Log( J R c ) 



V E2 + W) 



(6) 



The factor of 2 comes from the two parts of a symmetric geodesic, the 
— Log(i? c ) is a regulator, and r m is chosen to satisfy E 2 + f(r m ) = 0, so it 
is interpreted as the minimum of the r coordinate along the geodesic curve. 
Since the perturbation g(r) is localized near the singularity and decreases 
away from it, the maximum of g(r) is at r m . But r m must itself be quite 
large (of order 1 with our dimensionless setup) because the geodesic was 
taken to be far from the singularity, so g(r m ) must be extremely small. This 
gives a bound for L (from now on the regulator will be implicit) 

<L<2[™ ,„, * (7) 



which gives 

<L<Logf - ■ ^ , - - ) (8) 




-E 2 + l + \g(r m )\J ~ ~ *\-E* + l-\g(r. 

from which we see that L doesn't change much for geodesies far from the 
singularity, because such geodesies are characterized by small E, and in that 
range there is no sensitivity to the tiny |^(r m )| terms. In other words, the 
change in L is ~ |g(r m )|. 



Exponential Decay Models 

Now we will study a perturbation that decays exponentially as a function of 
the r coordinate. More precisely, we take g{r) = \ 2 e~ ar , where we will be 
taking a ^> 1. In this case, even at the black hole horizon, where r — 1, the 
perturbation will be exponentially small. 
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We are interested in correlation functions. In the large mass limit, we can 
use the WKB approximation, which is based on our knowledge of 



dr , , 

(9) 



VE 2 + r 2 - 1 - X 2 e~ ar 

There is no hope of doing this integral explicitly. It can be accurately approx- 
imated if it is split into two integrals, but this would involve us in complex 
algebra that can be avoided by a more direct approach. First note that in 
the standard, unperturbed the geodesic approaches the singularity, 

we have 

r°° dr 

L = 2 -> 2 log(r m ) (10) 

as r m — *■ 0. The exact solution in this case [3] gives L pa t as t — > oo, which 
is what we are seeing in the equation above. However, with the perturbation 
turned on, we get an integrand of the form (Ar 2 + Br + C)~2 from the 
expansion of the exponential at very small r, so both L and t are finite as 
r m -> 0. 

It is possible that we can recover a behavior similar to the unperturbed 
case by continuing r m past 0, and in fact, this is exactly what happens. The 
idea is this: we want to find an r m that takes the integrals for L and t to 
oo. Near this r m the integrand must behave like -^=. Thus we need to find 
a quadratic vanishing point of the function under the square root, namely 
E 2 + r 2 — 1 — A 2 e~ ar . This is easy to do, the result is that in the regime where 

aX < y|, we definitely have an r m value on the negative real r-axis. This is 
what we would like, since it simply means the true singularity has 'moved'. 
However, outside of this regime we still have an r m that allows for quadratic 
vanishing under the square root, but it will in general be complex, giving a 
complex E and generically complex L and t. Thus we will still find L pa ct 
for some c, but the physical interpretation is less clear. 
In the regime aA < 1 we can approximate using 

= E 2 + rl-l - X 2 e- ar - Pa A^r m + + C - ^ (11) 



by expanding the exponential to quadratic order and completing the square. 
Now if we set C = jj by tuning E we find that r m = — Jr. This allows us 
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to approximate the quadratic vanishing as 



L = 2 



Mr (12) 



+ Pe-) V / A(r + ^) 2 + C-g 



where R t is some small value of r where higher order terms from the ex- 
ponential become important. Both of these integrals do indeed diverge as 
we approach the r m noted above, and we can roughly approximate L pz 

^1 + 3 °g A4 ^ t for very large t by setting the rest of the integrand equal to its 

value at r m . Thus we get a linear L(t) for large t that is larger than in the 
unperturbed case. 

The next example is g(r) = X 2 e~ ar . We can use a similar approximation 
scheme, but we will get different results - L — > oo as r m — > 0, so the analysis 
is simplified. This is because when we expand the exponential we find 

(1 + X 2 a)r 2 m + (E 2 -l-X 2 ) = Ar 2 m + B^0 (13) 

so the perturbation does not destroy the quadratic vanishing at r m = 0. 
Again, it is possible to approximate the integral by breaking it up, but a 
simpler analysis gives the large t behavior. The important point here is that 
this behavior comes from the region very close to r m = 0, so we can ignore 
the rest of the integral. We find 

" fit Edr f Ri VI + \ 2 dr 



t = 2 



VAr 2 + B(r 2 -l-X 2 e- ar2 ) J Tm VAr 2 + B(l + A 2 - r 2 (l + X 2 a)) 

L (14) 



So we see that for very small r m corresponding to very large t, L{t) pa 
ty/l + A 2 , which is again larger than in the unperturbed case. 

A final question about the exponential decay models is at what value of 
t does the perturbed L(t) differ significantly from the original L(t)7 A quick 
way to get a rough estimate is to use equation (5) to compare an infinitesimal 
perturbed geodesic length element, 5s, to an unperturbed length element, 
5sq, both at a point r: 

5s r 2 — 1 — g(r) 



5s r 2 — 1 



(15) 
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Thus for a perturbation g(r) = \ 2 e~ arb , this expression is significantly dif- 
ferent from unity for r < a -1 / 6 . Here we interpret r as r m , the r of closest 
approach, so we see that the presence of g{r) starts to make a difference 
for = 1 — E 2 cT 2 l h . At this level of approximation we can use the 
unperturbed expression for t(E), giving 
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t, si Log(4afc - 1) (16) 

where t s is the value of t where the perturbed and unperturbed L(t) functions 
begin to separate. Even for very large a this is quite small, indicating that 
information about the geometry near the singularity is not difficult to recover 
from boundary correlators. 

We can also approach the problem by estimating AL, the difference be- 
tween the perturbed and unperturbed geodesic length, and comparing it to 
the exact, unperturbed answer. To do this, note that 

a r/^x „ f m dr f°° dr dr 

AL(E) = -2 / — : = +2 



ro 



VE 2 + r 2 - 1 J rm VE 2 + r 2 - 1 - X 2 e~ arb VE 2 + r 2 - 1 

(17) 

where r and r m are zeroes of E 2 + /(r) = for the unperturbed and per- 
turbed cases, respectively. 

We will bound the second term and then compute the first. A similar 
result will be important in the D = 5 case for different purposes. The second 
term is smaller than 



,°° dr dr nT I X 2 e~ 

2 / - = 2Log 1 + 



VE 2 + r 2 - 1 - \ 2 e- ar ™ VE 2 + r 2 - 1 
whereas the first term is just 



(18) 



l2H±^ j (19) 

Since they contribute with opposite signs and the absolute value of the first 
term is bigger, we see that (taking absolute values) 

2L °g ( m Tx2 -ar> 7 < ^L(E) < 2Log | (20) 
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and both bounds are reasonable approximations for AL itself. Now we are 

nplification 

log ( rM ;; rbm/2 ) 



interested in or upon simplification 



log(2r ) 



(21) 



Thus we see that for r m ~ a~~5 and so t ^ log(o), the perturbed L differs 
from the unperturbed by a sizable percentage, confirming our simpler method 
above. It is essential to remember that this method and the previous are only 
estimates. We have been calculating L as a function of r m , not t, so we can 
only get a rough idea of when the perturbation becomes important. 

Since both methods give an imprecise, but very encouraging estimate 
for the t value where the perturbation becomes important, we will also give 
an upper bound above which the perturbation must affect L significantly. 
We found the exact asymptotic behavior for Lit), and for both exponential 
models it was L ps ct for large t, for some c, where c = 1 in the unperturbed 
case. Once in this asymptotic regime, the perturbation is certainly impor- 
tant. This regime is characterized by an L dominated by the segment of the 
geodesic close to the singularity. Noting that for t = we have r m = 1 and 
L = 21og(2), we see that if L 3> 1, L must be dominated by contributions 
from inside the horizon. Since L oc t, for any t ^> 1 the perturbation makes 
a large effect, so in particular, for t m a the perturbation must be significant. 
Thus our previous analyses show that the perturbation should be important 
for t ~ log(a), and this calculation demonstrates that its effect will certainly 
have been long established for t m a. 

An Exactly Solvable Model 

The last perturbation to be considered for D = 3, g{r) = ^, is really a 
very large perturbation near the singularity, but it has the advantage that 
all integrals can be done exactly. Now we have 

f°° dr 
L = 2 (22) 

Jr m s J E 2 + r 2_ l _X 

The integral is doable though nontrivial, defining r m as root of the denomi- 
nator of the integrand, we get the final result 
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Note that as r m — > 0, the term implies that E — ► oo, implying that 
L — ► — oo. This is dramatically different from the unperturbed case because 
of the minus sign, but this is not surprising since the effect is huge near r = 0. 
The effect here is actually essentially the same as in the standard D = 5 
model, which is obvious given the similar behavior of this perturbation and 
the D = 5 model near the singularity. To see this more concretely, note that 

f°° Edr , . 

t = 2 (24) 

Jr m ( r 2_l_^)^2 + r 2_ 1 _ A 

As noted above, the limit r m — > corresponds to E — > oo, so in this limit 
the E dependent terms in the numerator and denominator of the integrand 
cancel, and we have 

f°° dr 

t=2 l (25) 

which is finite since there is no divergence from r near or oo. This gives 
a t c similar to that found by [5. in the D = 5 case, where there is a pole 
in the WKB correlator e~ mL and the behavior of the L(t) function changes 
qualitatively. We can evaluate t c explicitly, simplification gives 



2tt 



V2 + 8A 



+ VT+JX + y / v / TT4A - i^ (26) 



Model Stringy Effects in D = 5 

Now we will consider a different sort of perturbation in the D = 5 case. The 
goal here is to examine how boundary correlators change in a model where 
the curvature singularity at r = is 'resolved', so that curvature scalars are 
large and finite where they used to be singular. This could be considered as 
a model for stringy effects that smear the singularity, but it is not intended 
as an actual guess at what these effects are. Thus we will take 

f(r) = r 2 - (27) 

so that the e = case corresponds to standard AdS-Schwarzchild. In the 
standard case there exists a special value t c — | that corresponds to the 
geodesic becoming nearly null, so that L — > — oo and the two point correlator 
blows up. We will find that this t c value is 'resolved' when e 7^ 0, so that the 
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pole in the boundary correlator at t c splits into two poles, both off the real 
t axis. 

It will be relevant in the following to note that when t and L are re- 
garded as functions of E, the t c point corresponds to E — > oo. The integral 
representation of t(E) for the finite e case can be parameterized as 

r EdxV^r^ 

t= - ; 28 

Jr* (x - -^)Jx 3 + (E 2 + e 2 )x 2 + (E 2 e 2 - l)x 



x+e 



Note that as E — > -, the coefficient of x under the square root in the de- 
nominator vanishes. Since in this same limit we have r m = 0, the integrand 
blows up as i, so t — > oo, which never occurs in the unperturbed case. Thus 
for finite but very small e, t(E) increases steadily for a while, levels off near 
t c , then blows up to infinity for E near -. Note that this divergence comes 
from the behavior of the integrand near x — 0, but the only difference be- 
tween this integrand and the integrand for the L(E) function is the factor of 

(x ^-j) in the denominator and the factor of E in the numerator. These 

factors combine to give L(t) ~ - for E near j, or t — > oo. 

However, this leaves us with a mystery, because in the original model 
L(t) w 2Log(t — t c ) near t c . Where has this t c singularity gone? To answer 
this question we look at 

L(E) = (29) 
Jrl y/x 3 + (E 2 + e 2 )x 2 + (E 2 e 2 - l)x 

because we must have L(t) — > — oo at the new value(s) of t c . For very large 
E (much greater than i) the integral expression for L gets large, because 
there is a large range of x where the integrand is only decaying like x~ 1//2 . 

To see this explicitly, we can approximate L(E) very roughly in this 
regime as 

r 2 edx r°° dx 



v/^x 2 + (E 2 e 2 - l)x J e 2 Vx 2 + E 2 x + E 2 e 2 - 1 

(the lower limit of integration is — e 2 , not 0, because — e 2 is the largest root 
of E 2 + f(r) = 0) which gives 

L = ^Log(2i+2i£e)+^ r Log(2£e+2 v / £ 2 e 2 - l)-Log f2e 2 + E 2 + 2 v / 2~Pe 2 " 
EE V 

(31) 
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The first two terms go to for large E, but the third term goes to — oo. This 
divergence arises from the large x behavior of the integrand, which explains 
why t is finite for large E - its integrand includes an extra factor of -. Thus 
despite the singularity in t at E — j, it still makes sense to consider the 
limit E — > oo in our search for the missing pole, because for very large E, 
L diverges while t stays finite (close to the original t c but slightly complex), 
imitating the t c behavior of the unperturbed model. Calculating the full 
approximate value of t(E) is a bit lengthy, but the imaginary part of t c is 
easier because it must come from the region where x is negative. Thus we 
can get it from 



Im(£ c ) rs Im 




-Ee(x + e 2 ) 



2 y/E^x 2 + [E 2 t 2 - 1) 



x 



me 



3 



= ±— < 32 > 

where the ± comes from the choice of sign in the square roots of negative 
quantities. This shows how a perturbation that 'resolves' the black hole 
singularity splits the Log(t — t c ) pole. Without the perturbation we had 
L(t) ~ 2Log(t — t c ) near t c , so this must be modified at finite e. A good 
guess would be that L(t) Log((t — t c+ )(t — i c -)), but this remains an 
estimate because the integrals cannot be done algebraically. 



Perturbations and Power Series 

In this section we study f(r) = r 2 — 1+e r2 ar in the D = 5 case. In particular, 
we will analyze the power series for L(t) and show that the coefficients of 
this power series differ from the coefficients of the unperturbed case by an 
amount larger than the original coefficients after only a terms. This suggests 
that from a computational point of view our calculations are surprisingly sen- 
sitive to the geometry near the singularity. Though we choose this model for 
variety, our methods carry over straightforwardly to the D = 3 exponential 
models studied earlier. 

First we note that for small E, t(E) m E, so we will only worry about the 
L(E) function - studying small E is equivalent to studying small t. This is an 
approximation, because we are assuming that higher powers of t in E(t) do 
not conspire to destroy the results we find. First consider the unperturbed 



12 



L(E) function 



U E) = 2 I " , dT = Log | 1 (33) 

E 2 + r 2 -\ 




where the denominator of the integrand is zero at r . It will be useful later 
to note that the coefficient of the An th term in the expansion of r (E) is given 
approximately by Even more roughly, the coefficient of the n th power 
of E in the expansion of Lq(E) is about 2~ n / 2 . Now define 

f°° dr r°° dr 

AL(E) = 2 dr -2 (34) 

Jr m ^ E 2 +r 2_ 1±^I Jr ^£2 + r 2 _ £ 

where r m satisfies E 2 + f(r m ) = in the perturbed case. We will now examine 
the power series for r m ; this will be important because later we will show that 
the coefficients of the power series for AL(E) are essentially proportional to 
the coefficients in the series for r m (E). We can write 

r m = 1 + x + c k Ek ( 35 ) 
which we can plug into the defining equation to get 

E 2 (l + x + J2 c * Ek f + (! + x + J2 CkE ^ ~ 1 " e' a ~ ax e a ^ CkEk = (36) 

We can then use this equation to solve for the Ck by equating powers of 
E. Assuming x has been correctly computed the coefficients of all the odd 
powers of E will vanish (since the equation defining r m only involves E 2 ). 
We compute c 2 as 

E 2 + 4(1 + 3x)c 2 E 2 - e- a - ax ac 2 E 2 = 

c 2 w \ (37) 

Now we only need to keep the c 2 part when we expand the exponential, 
because terms involving higher c n will have fewer factors of a, which is a 
large number. By expanding the exponential we find 

e -a a n e -a a n 
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What we are interested in is the ratio of these terms to the terms in the 
expansion of L . This ratio is e n °°" , and it is maximized at the a th term, 
which is of order 1. 

Now looking back to equation (34) and noting that r m > r , we can write 

-dr f°° dr dr 



AL(E) = 2 / - aT +2 / 

Jro jE 2 +r 2 -\ Jr„ 



( / (39) 

We will now show that the contribution of the second term to the power series 
is exponentially suppressed as e~ a / 2 , and that the first term contributes to 
the power series of AL(E) terms proportional to (or larger than) the r m 
series. First we bound the second term. There is no theorem here that says 
that our bound on the second term also bounds the power series in E, but 
it should be extremely plausible in this case, because our bound does not 
destroy any violent dependence on E. In any case, we can always just regard 
our bound as an approximation. It is 

„a f°° dr dr 

2 nd term < 2 



'■'» yjE* + r*-l±^ ^E 2 + r 2 -± 



( p -ar m -2ar m \ 

where B = r m e~ aTrn ^ 2 and A = E 2 + 2^ + 213. Thus every term is suppressed 
by e _arm//2 , which proves the claim since r m is nearly 1 for t near 0. Now 
we need only show that the first term gives the much larger contribution 
claimed. This is not very hard, because we can just evaluate the first term 
directly, it is 

-* _ Log /g + 2ri + y r? +g , f „-A 



Jro 



E 2 + r 2 -\ V E2 + 2r o 



2r 2 V2r m - 2r 4r 2 (r m -r ) 

L0 & 1 / i 1 T~ 4 r ■■ 

V r o + r 5 1 + n 

which when expanded gives a power series in E with terms at least as large 
as the r m series, because r ~ 1. Thus we have shown that after about a 
terms the AL(E) power series has coefficients that are large compared to 
the Lq{E) coefficients. This is noteworthy because the perturbation itself is 
suppressed exponentially in a. 
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Path Integral Treatment 

All previous results have been obtained from the WKB approximation, which 
applies in the large mass limit. Thus it would be interesting to see how large 
effects emerge from a small perturbation in the full quantum mechanical 
treatment of the problem. Furthermore, in [Jj it was shown that there exists 
a dual description of the D = 3 AdS-Schwarzschild physics, depending on 
whether Lorentzian correlators are obtained from the Euclidean correlator 
through analytic continuation with respect to the Schwarzschild time coor- 
dinate or Kruskal time. Thus it should be possible to see the large effects in 
an explicitly outside the horizon calculation. We will see how this is done, 
and the key will be analytic continuation. 

We wish to calculate the two point function for a massive free scalar field 
in the perturbed S-AdS background. Specifically, we want to evaluate this 
two point function between field operators on opposite sides of the Penrose 
diagram, with the black hole between them, so that we can compare the 
effect of the perturbations with the semi-classical treatment already given. 
In order to do this, we will compute in the Euclidean spacetime, and then 
analytically continue in the t coordinate. The primary object we will consider 
is the Euclidean path integral 



I = / £>$ e 



" J 2 



^+(^-l- fl (r))(ft$) 2 + ^W* ! 



(42) 



Now for g(r) = A 2 e ar we can expand the argument of the exponential to 
first order in A 2 , to get 



(43) 

Note that the integrals are only evaluated outside the horizon, since the 
behind the horizon region doesn't even exist in the euclidean case. Now we 
have approximated the first order effect of the perturbation as an insertion 
of an operator smeared over spacetime. Thus the two point function can be 
written as 

<*(xi)*C«*)> - /$ (o:i)$o(^)e-^ (Aae -^[^ ( ^ +(W ]' 



2 



2*2 



$ (xi)$o(x 2 ) f 1 - J d 3 x 



AV 



_( r 2 _ iy. 
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The first term in the last line is obviously just the original, unperturbed two 
point function, and the second term is the first order effect of the perturba- 
tion. This is equivalent to a position dependent mass for the field, since it 
appears as the integral of a function multiplied with the first derivative of 
the field squared. It can be calculated in the same way standard Feynman 
diagrams are, meaning that we will 'contract' each exterior field with one 
of the squared differentiated fields. Just as in a standard field field theory 
setting, we use D(x, y), the unperturbed propagator, to contract fields. This 
gives the second term as 



A($(x 1 )$(x 2 )) w / d 3 x 



AV 



1 



1) 



;d t D(xi, x)d t D(x2, x) + d r D(xi, x)d r D(x 2 , x) 

(45) 



Now the D = 3 AdS-Schwarzschild propagator is given by an infinite sum, 
which arises from an infinite set of 'images' of the AdS propagator under the 
discrete identification of points that creates the singularity. To avoid this 
complication, we will be forced to work in the r + — > oo limit, in which the 
propagator takes the much simpler form 



D 



— 5— -cos(r + At) H cosh(r + A0) 



-2/14 



(46) 



For definiteness and to avoid long expressions, we will only deal with the 
d r D(xi,x)d r D(x2,x) term (the other term gives a result of the same order), 
and we will make a few simplifying approximations along the way. This term 
is given by 



2?r 



(It 



a 



dr2hl\ 2 ( 



cosh(0) - ^|=cosh(At 1 ) 



\ 



^ ^cosh((/>) 



cosh(Ati)) 



2/1++1 



:...au 



) 

(47) 

The contour C is from — T to —2ni — T, but when we analytically continue 
we will take T — > 00 and C will move to avoid the singularities. We can 



approximate r m as 1 + 



and Atj 



t-U 



Our goal here is to make it clear how the geodesic approximation gives 
a large contribution. For this purpose it suffices to derive an equivalent 
of the geodesic approximation from this operator calculation. The 0-angle 
dependence will not affect any of the interesting dynamics, so for simplicity 
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we set = 0. The key step now is to write the integrand as an exponential, 
giving 

[ dt [ dr Ah 2 ^ 2 e - m [ 21o s(('"-^'* 2 - lcosh (*-*i))('*-v / '* 2 - lcosh (*-*2)))] 

J C J r m 



x e 



ar 2 -log 1- I cosh(t-ti) 1 ^=cosh(t-t 2 ) 



(48) 



In the m — > oo limit the the second exponential is irrelevant compared to the 
first. Thus the integral should be approximated by looking at the stationary 
point of this first exponential. But the term being exponentiated is just 

log((r- Vr 2 - lcosh(t-ti))(r- Vr 2 - 1 cosh(t - 1 2 ))) w + L x _ 2 (49) 

where we define Lj_ x to be the geodesic length between the boundary point 
i and the integration point x. One easy was of seeing this is that 

' e -mlog(r-Vr 2 -lcosh(At)) ^ g-mL (50) 



(r- v / r^Tcosh(At)) 2?i + 

Therefore taking the stationary point of this term is entirely equivalent to 
taking the point x to lie along the unperturbed geodesic connecting the two 
boundary points. We conclude that in the large m limit we should inte- 
grate the perturbation along the geodesic. This gives the same result as the 
geodesic calculation to first order in A. Since the perturbed two point func- 
tion is analytic in At, we see that the full quantum calculation must give the 
same results as the geodesic computation, including the large effect of the 
perturbation for small t. 



Conclusions 

We have seen in a variety of situations that data about the geometry near 
the spacelike singularity of the AdS-Schwarzschild spacetime is contained in 
boundary correlators. In each case, the recurring theme was that this in- 
formation can be calculated more easily than might be expected, and that 
it is sensitive to specifics such as the function perturbing or 'resolving' the 
singularity. This information is rather indirect, since we still do not under- 
stand the spacelike singularities themselves. However, it is possible that our 
method could be a clue to a more complete description, and it is exciting to 
see explicit signatures of the singularity. 
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